Simulation of Transitions between "Pasta" Phases in Dense Matter 



in 
o 
o 

(N 
00 

m 
> 

o 

00 

o 
o 



X 



Gentaro Watanabe"'^'"^, Toshiki Maruyama'', Katsuhiko Sato''^'^, Kenji Yasuoka-'' and Toshikazu Ebisuzaki^ 
"NORDITA, Blegdamsvej 11, DK-2100 Copenhagen 0, Denmark 
''Department of Physics, University of Tokyo, Tokyo 113-0033, Japan 
'^The Institute of Chemical and Physical Research (RIKEN), Saitama 351-0198, Japan 

''ASRC, Japan Atomic Energy Research Institute, Tokai, Ibarakt 319-1195, Japan 
"Research Center for the Early Universe, University of Tokyo, Tokyo 113-0033, Japan 
^Department of Mechanical Engineering, Keio University, Yokohama, 223-8522, Japan 

(Dated: February 8, 2008) 

Calculations of equilibrium properties of dense matter predict that at subnuclear densities nuclei 
can be rodlike or slablike. To investigate whether transitions between phases with non-spherical 
nuclei can occur during the collapse of a star, we perform quantum molecular dynamic simulations 
of the compression of dense matter. We have succeeded in simulating the transitions between rodlike 
and slablike nuclei and between slablike nuclei and cylindrical bubbles. Our results strongly suggest 
that non-spherical nuclei can be formed in the inner cores of collapsing stars. 
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In ordinary matter, atomic nuclei are roughly spherical 
because, in the liquid drop picture of the nucleus, effects 
of the nuclear surface tension are greater than those of the 
Coulomb forces. When the density of matter approaches 
that of atomic nuclei, calculations predict that, in equilib- 
rium state, the nuclei will adopt different shapes, such as 
cylinders and slabs, etc. These phases with non-spherical 
nuclei are often referred to as "pasta" phases 

In the initial stage of the supernova explosions, mat- 
ter in the collapsing iron cores experiences an adiabatic 
compression, which leads to an increase of the density in 
the central region from ^ \{f g cm^"^ to around the nor- 
mal nuclear density po = 0.165 fm^'^ just before the star 
rebounds; the temperature there reaches the order of 1 
MeV. The pasta phases are thus expected to be formed 
in the inner cores during the collapse of stars. How- 
ever, such a speculation is based on phase diagrams of 
the equilibrium state (e.g., Refs. 0,0 for finite temper- 
atures) or static and perturbative calculations 0, 0] . It 
is still unclear whether or not the pasta phases can be 
formed and the transitions between them can be realized 
during the collapse, which lasts less than a second. Be- 
cause of the drastic changes of nuclear shape that occur 
under non-equilibrium conditions, this problem is more 
difficult than the realization of the pasta phases by cool- 
ing at constant density, as demonstrated in Ref. [2| , and 
an ab-initio approach is called for. 

In the present Letter, we solve the problem about 
the transitions between pasta phases using a dynami- 
cal framework for nucleon many- body systems called the 
quantum molecular dynamics (QMD) [Sj . QMD is a suit- 
able approach to describe thermal fluctuations and is ef- 
ficient enough to treat large systems consisting of sev- 
eral nuclei. Furthermore, at the relevant temperatures 
of several MeV, shell effects, which cannot be described 
by QMD, are less important because they washed out by 
thermal fluctuations. 

The pasta phases have recently begun to attract the 



attention of researchers (see, e.g., Ref. [2j and references 
therein). The mechanism of the collapse-driven super- 
nova explosion has been a central mystery in astrophysics 
for almost half a century. Previous studies suggest that 
the revival of the shock wave by neutrino heating is a cru- 
cial process. As has been pointed out in Refs. 7, .1,QJ and 
elaborated in Refs. [ill ll^ . the existence of the pasta 
phases instead of uniform nuclear matter increases the 
neutrino opacity of matter in the inner core significantly 
^ due to the neutrino coherent scattering by nuclei 
3 ^1 ; this affects the total energy transferred to the 
shocked matter. Thus the pasta phases could play an im- 
portant role in the future study of supernova explosions. 

In the present study, we use a nuclear force given by a 
QMD Hamiltonian with medium-EOS parameter set in 
Ref. This Hamiltonian contains the momentum de- 
pendent "Pauli potential", which reproduces the effects 
of the Pauli principle phenomenologically. Parameters 
in the other terms of the Hamiltonian are determined to 
reproduce the saturation properties and the properties 
of finite nuclei in the ground state, especially of heav- 
ier ones relevant to the present study It is also 
confirmed that a QMD Hamiltonian close to the present 
model provides a good description of nuclear reactions in- 
cluding the low energy region (several MeV per nucleon) 
[IT'I, which would be important for the present case. 

Using the above QMD Hamiltonian, we perform simu- 
lations of symmetric nuclear matter with 16384 nucleons 
in a cubic box with periodic boundary condition (see Ref. 
12] for other cases). The system contains equal numbers 
of protons (and neutrons) with spin up and spin down. 
The relativistic degenerate electrons which ensure charge 
neutrality are treated as a uniform background because, 
at subnuclear densities, the effect of the electron screen- 
ing is small . The Coulomb energy, taking account 
of the Gaussian charge distribution of the proton wave 
packets, is calculated by the Ewald method. The temper- 
ature T is measured by the effective kinetic temperature 
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for momentum-dependent potentials, which is consistent 
with the temperature in the Boltzmann statistics ■ The 
QMD equations of motion are integrated by the fourth- 
order Gear predictor-corrector method with a multiple 
time step algorithm. Integration time steps At are adap- 
tive in the range of At < 0.1 - 0.2 fm/c. 

As the initial condition, we use samples of the colum- 
nar phase and of the laminar phase of 16384-nucleon 
system at T ~ 1 MeV obtained in Ref. 0|. In prepar- 
ing them, we first combine eight replicated 2048-nuclcon 
samples in the ground state (T ^ MeV; nucleon number 
density p = 0.225po for the phase with rodlike nuclei and 
0.4po for the case of slablike nuclei) into a 16384-nucleon 
sample, and then put random noise on the positions and 
the momenta of nucleons up to 0.1 fm and 1 MeV/c, re- 
spectively. We equilibrate the sample at T = 1 MeV for 
^ 4000 - 5000 fm/c using the Nose-Hoover thermostat 
for momentum-dependent potentials 0. We further re- 
lax the sample for 5000 fm/c without the thermostat. 

Starting from the above sample, we simulate the adi- 
abatic compression. In the case starting from the phase 
with rodlike nuclei [slablike nuclei], the density is in- 
creased by 2xl0~''po [lxlO~*/Oo] every 100 steps by 
changing the box size (the particle positions are rescaled 
at the same time). The average rate of change for the 
density is ~1.3x 10-Vo/(fm/c) [-7.1 x 10-Vo/(fm/c)]; 
this rate ensures the adiabaticity of the simulated com- 
pression process with respect to the change of nuclear 
structure Finally, we relax the compressed sample 
at p =0.405 po [0.490 po]. These final densities are those 
of the phase with slablike nuclei [cylindrical bubbles] in 
the equilibrium phase diagram at T ~ 1 MeV jj] . 

The resulting time evolution of the nucleon distribu- 
tion is shown in Figs. 1 and 2. Starting from the phase 
with rodhke nuclei [Fig. 1-(1); p ~ 0.225po (volume frac- 
tion of nuclear matter region u — 0.20 - 0.22 at time 
t = 0], the phase with one-dimensional layered lattice of 
slablike nuclei is formed [Fig. l-(9); p = 0.405po {u = 0.41 
- 0.45) at < = 17720 fm/c]. During the compression, 
the temperature increases gradually up to ~ 1.35 MeV 
in the final state. We note that, in the process of the 
compression, the phase with rodlike nuclei persists as 
a metastable state, and moreover, until nuclei begin to 
touch and fuse they are not elongated along the plane of 
the final slabs [see Fig. l-(2)]. This shows that the tran- 
sition from the phase with rodlike nuclei to the slablike 
nuclei is not triggered by the fission instability. This re- 
sult is consistent with a previous study, which shows the 
stability of the rodlike nuclei against a small quadrupo- 
lar deformation of the cross section When the in- 
ternuclear spacing becomes small enough and once some 
pair of neighboring rodlike nuclei touch due to thermal 
fluctuations, they fuse [see the lower two nuclei in the 
middle column in Figs. l-(3) and l-(4); u = 0.27 - 0.30 
at t = 6050 fm/c]. Like a chain reaction, such connected 
pairs of rodlike nuclei further touch and fuse with neigh- 



boring nuclei in the same lattice plane [see Figs. l-(5) and 
l-(6)]. Each fusion process in the chain reaction proceeds 
on a time scale of order 10^ fm/c, which is much shorter 
than the time scale of the density change . 

The transition from the phase with slablike nuclei 
to the phase with cylindrical holes is shown in Fig. 2 
(u = 0.42 - 0.45 and 0.55 - 0.59 at i = and 27370 
fm/c). When the internuclear spacing decreases enough, 
neighboring slablike nuclei touch due to the thermal fluc- 
tuation as in the above case. Once nuclei begin to touch 
(u = 0.49 - 0.52 at t = 8460 fm/c), bridges between the 
slabs are formed at many places on a time scale (of or- 
der 100 fm/c) much shorter than that of the compression 
[cf. Figs. 2-(3) and 2-(4)]. After that the bridges cross 
the slabs nearly orthogonally for a while, which makes 
hollow regions on a square lattice rather than the flnal 
triangular one. Nucleons in the slabs continuously flow 
into the bridges, which become wider and merge together 
to form cylindrical holes. Afterwards, the connecting re- 
gions consisting of the merged bridges move gradually, 
and the cylindrical holes relax to form a triangular lat- 
tice. The final temperature in this case is ~ 1.3 MeV. 

Let us now investigate the detailed time evolution of 
the nuclear structure. The integral mean curvature and 
the Euler characteristic (see, e.g., Ref. ^3]) powerful 
tools for this purpose. Suppose there is a set of regions R, 
where the density is higher than a threshold density pth- 
The integral mean curvature and the Euler characteristic 
for the surface of this region dR are defined as surface in- 
tegrals of the mean curvature H and the Gaussian curva- 
ture G, respectively; i.e., /^^ HdA and x = ^ Jg^j GdA, 
where dA is the area element of the surface of R. The 
topological quantity x = (number of isolated regions) — 
(number of tunnels) + (number of cavities) ^2^] . We cal- 
culate these quantities using the density field of nucleons 
on 128^ grid points (see Ref. for detailed procedures). 

In Figs. 3-(A) and 3-(B), the quantities Jgj^HdA and 
X calculated for each nucleon distribution in Fig. 1 are 
shown as functions of pth- In this case, the initial con- 
dition is the phase with rodlike nuclei, whose structure 
is well characterized by the plateau of J^^ HdA ~ 4000 
fm [see the red curve in Fig. 3- (A)]. The slablike nuclei 
in the final state, on the other hand, are characterized 
by Jgj^ HdA ~ fm, corresponding to the plateau of the 
gray curve in Fig. 3-(A). 

The behavior of x clearly shows that the rodlike nuclei 
begin to touch between t = 5290 and 6050 fm/c, when the 
plateau value of x starts to deviate from zero, character- 
izing the rodlike nuclei, to negative values, characterizing 
the multiply connected structures such as sponges. We 
note that before the nuclei touch, the change in /g^ HdA 
is small except for lower values of pth ^ O.lpo- This re- 
flects the facts that the phase with rodlike nuclei persists 
as a metastable state and that the transition is not in- 
duced by the flssion instability Also the behavior 
of X should be noted; as can be seen from Fig. 4-(A), it 
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FIG. 1: (Color) Snapshots of the transition process from the 
phase with rodhke nuclei to the phase with slabhke nuclei 
(the whole simulation box is shown). The red particles show 
protons and the green ones neutrons. After neighboring nuclei 
touch as shown by the circle in Fig. l-(3), the "compound 
nucleus" elongates along the arrow in Fig. l-(4). The box 
size is rescaled to be equal in this figure. 
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FIG. 2: (Color) The same as Fig. 1 for the transition from the 
phase with slablike nuclei to the phase with cylindrical holes 
(the box size is not rescaled in this figure). After the slablike 
nuclei begin to touch [see the circle in Fig. 2-(3)], the bridges 
first crosses them almost orthogonally as shown by the arrows 
in Fig. 2-(5). Then the cylindrical holes are formed and they 
relax into a triangular lattice, as shown by the arrows in Fig. 
2-(9). 
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FIG. 3: (Color) The quantities J^^ HdA and x functions 
of the threshold density pth calculated for the nucleon density 
fields of Fig. 1 [(A) and (B)] and of Fig. 2 [(C) and (D)]. 
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FIG. 4: (Color) Time evolution of HdA and x during 
the simulations. The data points and the error bars show, 
respectively, the mean values and the standard deviations in 
the range of pth = 0.3 - 0.5po, which includes the values cor- 
responding to the nuclear surface assumed in (20t . We set 
the averaging range relatively wide because the nuclear sur- 
face cannot alw ays be characterized by a single value of pth 
as mentioned in '2(J|. The panel (A) is for the transition from 
cylindrical [C] to slablike nuclei [S] and the panel (B) for the 
transition from slablike nuclei to cylindrical holes [CH] . Tran- 
sient states are shown as (C,S) and (S,CH) for each transition. 



becomes negative between the phase with rodlike nuclei 
(x ~ and HdA > 0) and the phase with slabhke 
nuclei (% ~ and Jg^ HdA ~ 0). This implies that the 
transition proceeds through a transient state with "spon- 
gelike" structure. The state which gives the smallest x at 
t ~ 9840 fm/c corresponds to the moment when all of the 
rodlike nuclei are connected to others by small bridges; 
after that, the connected nuclear rods relax into slablike 
nuclei, i.e., the bridges in the slablike structures merge 
to form the nuclear slabs and those across the slabs dis- 
appear. The whole transition process can be divided into 
the "connecting stage" and the "relaxation stage" before 
and after this moment; the former starts when the nuclei 
begin to touch and it takes ~ 3000 - 4000 fm/c and the 
latter takes more than 8000 fm/c. 
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The same quantities are shown for the transition from 
the phase with slabhke nuclei to the phase with cyhndri- 
cal bubbles in Figs. 3-(C), 3-(D), and 4-(B). The initial 
and the final structures are characterized by the plateau 
values of J^^ HdA ~ and ~ —2000 fm, respectively. 
From Figs. 3-(D), and 4-(B), we see that the slablike nu- 
clei begin to touch at t 8460 fm/c and the connection 
of the slablike nuclei by the small bridges are completed 
at t < 12000 fm/c corresponding to the state with the 
lowest X- III this case, the connecting stage lasts for 
~ 3000 - 4000 fm/c and the relaxation stage for more 
than 15000 fm/c. In the latter period, the bridges merge 
to form cylindrical holes shown by the increase of x tO" 
ward zero, and, simultaneously, their positions relax into 
a triangular lattice as mentioned before. 

In conclusion, we have succeeded in simulating the dy- 
namical process of two types of transitions between pasta 
phases at sub nuclear densities. Our calculations support 
the idea that transitions between pasta phases can oc- 
cur during stellar collapse. The particular transitions we 
have examined are triggered by the thermal fluctuation, 
not by the fission instability. They consist of the con- 
necting stage and the relaxation stage. The total time 
of the connecting stage is 3000 - 4000 fm/c in our sim- 
ulations, which could be shortened by the artificial com- 
pression. However, we can conclude that the connecting 
stage would be complete in a time scale of order 10'^ fm/c 
taking account of the facts that each connecting process 
observed in this stage proceeds much faster than the com- 
pression and that the time scale of ordinary nuclear fis- 
sion is about 1000 fm/c. The relaxation stage takes about 
10000 fm/c or more. A remaining challenge is to investi- 
gate the transition from the bcc lattice of spherical nuclei 
to the triangular lattice of rodlike nuclei 12^. If this pro- 
cess is confirmed, the existence of the pasta phases in 
supernova cores will be almost established. 
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